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We develop a transfer matrix method to compute exactly the spin-spin correlation 
functions < s^Sn > of Bethe lattice spin models in the external magnetic field h and 
for any temperature T. We first compute < SQSn > for the most general spin - S Ising 
model, which contains all possible single-ion and nearest-neighbor pair interactions. 
This general spin - S Ising model includes the spin-1/2 simple Ising model and the 
Blume-Emery-Griffiths (BEG) model as special cases. From the spin-spin correlation 
functions, we obtain functions of correlation length, ^(T, h), for the simple Ising model 
and BEG model, which show interesting scaling and divergent behavior as /i ^ 
and T approaches the critical temperature Tc- Our method to compute exact spin- 
spin correlation functions may be applied to other Ising-type models on Bethe and 
Bethe- like lattices. 



PACS numbers: 05.50+q, 75.10-b 

I. INTRODUCTION 

The Bethe |l| and the Bethe-hke lattices have been widely used in solid state and 
statistical physics P-[TE[| as they represent underlying lattices for which many problems can 



be solved exactly. The Bethe lattice has attracted particular interest because it usually 
reflects essential features of systems, even when conventional mean - field theories fail fill . 



Besides the thermodynamic quantities, such as the magnetization, the specific heat, etc., 
the correlation function contains important information about a phase transition system 
T9| and is often studied by theoretical calculations pO|-p5|] and experimental measurements 



^ . It is widely believed that the singular behavior of physical quantities at the critical 
temperature Tc of a second order phase transitions is related to the divergence of the corre- 



lation length ^ at Tc [|T9| . In connection with this, the knowledge of the exact form of the 
spin-spin correlation function is very crucial for locating phase transitions and for analytical 
investigations of physical phenomenon. Furthermore, for many experiments, spin-spin corre- 
lation functions are most relevant, since they are measured by standard probes such as linear 
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response to an adiabatic or isothermal applied field, or scattering of neutrons or electromag- 
netic waves |26|. In the past several years, considerable progress has been achieved in the 
computation of the correlation function of statistical mechanical systems p0|,p2|. However, 



the spin-spin correlation functions are exactly known only in a few models, including the 
two-dimensional Ising model in zero magnetic field and at any temperatures |21| , the planar 
Ising model in an magnetic field and exactly at Tc 



22 



A long-standing problem of statistical mechanics is the exact solution of the spin-spin 
correlation function for the Ising model in an external magnetic field and at any tempera- 
ture. In this paper, we develop a transfer matrix method to compute exactly the spin-spin 
correlation functions < s„So > for the most general spin - 5* model on the Bethe lattice for 
any temperatures T and external field h. The model include the Ising model and the 
Blume- Emery-Griffiths (BEG) model |^ as special cases. The correlation length C,{T,h) 
obtained from the spin-spin correlation function shows interesting scaling and divergent be- 
havior as h ^ and t — 0, where t = (T — Tc)/Tc [^. Our results thus solve a long 
standing puzzle in the critical phenomena of the Bethe lattice Ising model. Our method 
may be applied easily to other Bethe lattice spin models. 

Although the free energy of the Cayley tree Ising model in zero external magnetic field is 
an analytic function of the temperature, the magnetization m and the magnetic susceptibility 
X of the central spin sq of the Bethe lattice Ising model have singular behaviors with the 
critical exponents (3=1/2, 6=3 and a=0 However, there is no previous calculation 



which shows that ^ of the Bethe lattice Ising model has singular behavior. Our work shows 
clearly that ^ diverges at the critical point of m and x fo'^ the Bethe lattice Ising model. 

The difference between the Cayley tree and the Bethe lattice has been discussed by 
Baxter |30]. In the Cayley tree, the surface plays a very important role because the sites 
on the surface is a finite fraction of the total sites even in the thermodynamic limit. As 
a consequences, the spin models on the Cayley tree exhibits quite unusual type of phase 
transition without long - range order pB|-p5|; the calculated correlation functions do not show 
singular behavior p3|-p5|]. To overcome this problem one usually considers only properties 
of sites deep in the interior (away from the surface) of the Cayley tree. The union of such 
equivalent sites, with the same coordination number q, can be regarded as forming the Bethe 
lattice |3^. Thus the Bethe lattice is assumed to have translation symmetry like any regular 
lattices. 

In this paper, we demonstrate the crucial role of Bethe lattice dimensionality in deter- 
mining critical behavior of the correlation length and show clearly that correlation length 
^ diverges at the critical point of m and x fo^' the Bethe lattice Ising model with critical 
exponent = 1, which is different from the mean-field critical exponent u = 1/2 considered 
by Tsallis and Magalhes in a recent review paper |jlO[, but is consistent with the critical 
exponent of the localization length associated with density-density correlator in the Bethe 
lattice Anderson model obtained by a supersymmetry method 0]. Our result gives indepen- 
dent support to the idea that mean - field approximation and Bethe lattice approach are not 
equivalent in principle [TT]. In this paper we also analyze exact spin - spin correlation func- 
tion for the BEG model ||3T| and find that at the tricritical point, ^ ~ with tricritical 
exponent ut = 1/2. 

Thus, we have solved a long standing puzzle in the critical phenomena of the Bethe lattice 
Ising model. Our approach can be extended easily to other models (Potts, Ashkin- Teller, 
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etc.) on the Bethe and Bethe-like lattices. In particular, our results for the Ising model on 
the Bethe lattice can be extended easily to Husimi lattices, because they can be related with 
each other in terms of the star - triangle transformation [pO| . 

The outline of the paper is as follows: In Sec. II the most general spin - S model is 
defined and the spin - spin correlation function of this model is evaluated in a closed form. 
In Sec. Ill we analyze the exact correlation function of the Bethe lattice spin - 1/2 Ising 
model and discuss the critical properties. In Sec. IV the critical behavior of the spin - spin 
correlation function is analyzed for the spin - 1 Ising model on the Bethe lattice. In Sec. V 
a brief discussion of our results is presented. 



II. A MOST GENERAL SPIN-5 MODEL 



The Ising - type models with spin greater than 1/2 have rich fixed-point structures. The 
great interest in these models arises partly from unusually rich phase transition behavior 
they display as their interaction parameters are varied, and partly from their many possible 
applications. 

A spin-1 Ising model, which was initially introduced by Blume-Emery-Griffiths (BEG) 
?T| in connection with phase separation and superfluid ordering in He^ - He^ mixtures. 



The BEG model has played an important role in the development of the theory of multi- 
critical phenomena associated with various physical systems |^3| and has extensively been 
investigated in the hterature 

The spin-3/2 Ising model with dipolar and quadrupolar interactions was first introduced 
to explain phase transitions in DyVO^ by Sivardiere and Blume and a different spin-3/2 



Ising model for ethanol- water-carbon-dioxide was introduced by Krinsky and Mukamel ||35|| . 
Another spin-3/2 Ising model was investigated by Sa Barreto and De Alcantra Bonfim [3U . 
Let us define the most general spin-S" model by the Hamiltonian 

-PH=Y: H^{s„s,) + J2H,{s,), (1) 

<ij> i 

where P = {ksT)'^ and Si is a spin variable which takes a value on {—S, —S+1, S' — 1, S}. 
The first sum goes over all nearest - neighbor pairs of the Bethe lattice and the second over 
all sites. Hi{si,Sj) contains all possible nearest - neighbor pair interactions and can be 
written as 

25 

Hi{s„s,)= ■J.us^^",, (2) 

fj,,u=l 

H2{si) includes all possible single ion interactions 

25 

HM) = E K^t- (3) 

There are S{2S + 1) independent nearest - neighbor coupling constants {J^u) in Eq. (0) and 
25* external fields in Eq. @. The Hamiltonian (|lD can describe a variety of different 
model. The partition function has the form 



3 



Z = J2exp\j2 H^{s.,s,) + Y.H2{s,)\. (4) 

{s} [<ij> i ) 

The advantage of the Bethe lattice is that for models formulated on it exact recursion 



relations can be derived. The calculation on a Bethe lattice is done recursively |30|. When 
the Bethe lattice is "cut" apart at the central site 0, it separates into q identical branches, 
each of which contains (g — 1) branches. Then the partition function of the model can be 
written as 

= J2exp{-PH) = 5:exp{/f2(3o)K(so), (5) 

s so 

where sq is a spin in the central site, N is the number of generations (A^ oo corresponds 
to the thermodynamic limit), and gN{so) is in fact the partition function of one branch. 
Each branch, in turn, can be cut along any site of the Ith-generation which is nearest to the 
central site. The expression for gN{so) can therefore be written in the form 

9Niso) = Y^exp {Hi{so, si) + H2{si)}gJ^_^{si), (6) 

where 7 = g — 1. 

Consequently, we have 25" + 1 recursion relations for the 5'7v(so), where sq takes values 
{—S, —S + 1, S — 1, S). After dividing each of recursion relations on a recursion relation 
for gN{S), we have 25" recursion relations for xn^Sq) 

, . ^ Esi exp {ij'i(3o, si) + ^2(^1)} /„v 

''^^'"^ Es,exp{iJi(^,si) + if2(si)}x],_i(si)' ^ ^ 

where 

2:iv(so) = fi'Jv(so)/fi'Af(5') (8) 

and the equation for gN^S) is 

9n{S) = gl_,{S)J2exp {H^{S, s,) + H2{s^)}xl_,{s,). (9) 

Since the right - hand side of Eq. (|^ is bounded in x^, it follows that x^ is finite for 
^ 00. Through x^, obtained by Eq.(^, one can express the density =< (^^^ > of 
central site (the symbol < ... > denotes the thermal average), where n take values from 1 
to 2S: 

' \SJ ^ Es,exp{H2{so)}x%{so) ' ^ 

and other thermodynamic parameters. 

So we can say that the x^ in the thermodynamic limit (A^ — > 00) determine the states 
of the system. For this reason the recursion relations for x^ given by Eq. (|^) can be 
called the equations of state for the spin - 5* model on the Bethe lattice. For example, at 



4 



high temperatures the recursion Eq.(0) tends to a fixed point and therefore the system has 
an appointed magnetization m. Equations (|^) and ([l^) are fundamental equations for the 
spin-S" model on the Bethe lattice. 

For calculating the spin - spin correlation function, it is now convenient to write down 
the expression for partition function in the following form 



Zn= Y1 exp ^ i^i(Si,Si+i) +X]^2(Si) ]9N{sQ)gN-l{Sl)---9N-n+l{Sn-l)gN-n{Sn), 



SoSl...S„ \j=0 2=0 



(11) 



where n denotes the number of steps from the central point 0. Summing over s„, s„_i, si 
consistently, we obtain Eq.(^ again. Then, the two-spin correlation function between sq 
and s„, r(?T,) = ^ < SQSn >, can be written as 

(12) 

We will show below that the calculation of < SqS^ > can be performed by a transfer - matrix 
method. Other techniques for the calculation of the spin - spin correlation function are not 
able to give the correlation in the presence of an external magnetic fields. For example. 



the results obtained in p3|-p5|] are appropriate only for h strictly equal to zero without any 
symmetry breaking effects. It should be noted that in order to obtain results which are 
relevant for the Bethe lattice we must used the proper thermodynamic limit N oo. 

The properties of the Bethe lattice can be investigateded by considering a Cayley tree 
with a very large number of generations (A^) and one only looks at the thermal ensemble 
of the sites in the interior part of the first n generations. Then the limit — >■ cxd is taken 
before n — > cxd. 

We are interested in the case when the series of solutions of recursion relations given by 
Eq. (0) converges to a stable point as ^ oo. In the thermodynamic limit (A^ oo), we 
may expect x^^n does not depend on n, so that XAr_„ and can be regarded as the same 
fixed - point solutions x of recursion relations given by Eq. (|^), which corresponds to the 
behavior in the interior part of an infinite Cayley tree (i.e. the Bethe lattice). In this case 

lim X]\f_n{s) = x{s) 

for all finite n. Then the recursion equations (equations of state) becomes 

. X _ Es^exp {Hi{so, Si) + H2{Si)} x^jsi) 

Es,exp{H,{S,s,) + H2{s^)}x^{s^) ^ ^ 

and the spin-spin correlation function T{n) given by Eq. (|T^ takes the form 

^^^^^ _ 1 EsQS...... ^^^^'•^'+^^+^-^^^^^'^)x^(^0)x^-^( (^4^ 
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The correlation function of Eq. ([M]) can be expressed in the vector-matrix form. For 
this purpose, let us introduce a (25" + l)x(2S + 1) matrix V and a (25* + l)-component 
column vector R. The elements of V are given by 



Vss' = exp \Hi{s, s') + 



H2(s) + H2(s') 



\x(s)x(s')]''^' 



(15) 



where s and s' independently take values S,S — 1,...,— 5 + 1,—S. The vector R and 
transposed vector R-^ have elements 



Ho(s) 



exp 



1x^2 (s). 



(16) 



Let us also introduce the diagonal matrix S 



/I ... \ 
^ ... 



VO ... -IJ 

With these definitions we may rewrite Eq. (|T^) in the vector-matrix form 

R^SV"SR 



r(n) 



R^V"R 



(17) 



The transfer matrix V is real-symmetric {Vss' = K's), and can be diagonalized by the 
transformation 

/Ai ... \ 
As ... 



P^VP = 

V . . . A25+1 J 

where P is a (25* + 1) x (25* + 1) matrix with the elements Pss' and Ai , A2, ... and A25+1 - are 
the eigenvalues of the matrix V. These eigenvalues can be obtained from the characteristic 
equation. Then the spin-spin correlation function of Eq. ([T7D can be written as 



Tin) 



R^SPLP iSR 



(18) 



R^PLP iR 

where L is a (25* + l)x(25' + 1) diagonal matrix 



fX'^ ... \ 
A^ ... 



V . . . A^^+i J 

After some algebraic manipulations, using Eqs. (p!3D, (p!5|) and (|T6|), we may write the 
correlation function of Eq.(|l^) as 
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2S 

T{n)=ml + Y.A,ll (19) 

k=l 

where 4 = Xk+i/Xi- An explicit expression for has been given in Appendix. It should be 
noted that the obtained exact expression for the spin - spin correlation function depends on 
the ratios of the eigenvalues Afc/Ai. Thus, Eq. (|19]) together with Eqs. (p!OD and (|1^) give 
us full set of equations for investigation spin-S* model on the Bethe lattice. In the following 
we turn to various examples. 

III. SPIN-1/2 ISING MODEL 

We first consider a spin-1/2 Ising Hamiltonian at a temperature T and an external 
magnetic field h 

- (3H = AJ SiSj + 2hY,Si, (20) 

<ij> i 

where Sj = +1/2 or —1/2 and first term describes the ferromagnetic coupling (4J) between 
the spin at site i and j. 

For the magnetization (mi = 2 < sq >) of the spin in the central site we can obtain from 
Eq. (|10D the following expression 

exp (2h) — x'^ . 
^1 = ) , i , (21) 

where x is the fixed point of the recursion relations ( pTSD in the thermodynamic limit 

.«-exp(-2ft) = i^^EM^. (22) 

exp (2 J) — X 

The two-spin correlation function for the spin- 1/2 Ising model can be obtained from Eq. 
(0) 

V{n) = m\^A^{^ , (23) 

with Ai = 1 — ml (see Appendix) and Ai, A2 - are the eigenvalues of the 2x2 matrix V. 
The elements of V are given by 

Vss' = exp {AJss' + hs + hs) \x{s)x{s')]^ , (24) 



where s and s' independently take values ±1 

V = 



The eigenvalues of the matrix V can be obtained from the characteristic equation 

A^ - A(y++ + y__) + -1/2 = 0. (25) 
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Using Eqs. (|2|) and (HD we find that 

Ai = 2 sinh (2 J) 



exp (J + h) 
(exp (2 J) — X 



Aa = (2 cosh (2 J) 



exp (J + h) 



x] 



(exp (2 J) 

Thus, the spin-spin correlation function can be expressed exactly as 

r{n) =m2 + (l-m2)A", 

where 



2 cosh (2 J) — X — x ^ 
2 sinh (2 J) 



(26) 



(27) 



(28) 



mi is the magnetization of the spin in the central site given by Eq. (|2lD and x is the solution 
of the recursion relation given by Eq.(|2^. 

It is well known that the Ising model on the Bethe lattice exhibits ferromagnetism, with 
a critical point at /i = 0, x = 1, T = T^, where Jc = | In and critical exponents (3, 6 and 
a have the "classical" values /? = |, 5 = 3 and a = 13^]. Let us now consider the general 
behavior of the spin-spin correlation function in the critical region. 

First consider the case h = and T = Tc. From Eq. (p7|) , we obtain T{n) = {q — 1)^". 
Let us consider a Bethe lattice with coordination number q, whose dimension is defined by 
dn = which tends to infinity with n ^ oo for g > 2 and equal 1 for q = 2, where 
r< — 

becomes the ordinary one- dimensional chain. In the limit of large n, dn for all g > 2 becomes 

n 



„ is the total number of sites. We should note that for g = 2 the Bethe lattice 
ordinary one- dimensional chain. 

ln(g-i: 



d 



Inn 



Thus we can write, for large n, 



T{n) = {q-1] 

Near the critical point, setting as usual t - 
function is 



r(n) 

where the correlation length ^ is given by 

1 



exp 



Tc 



n-"^. (29) 
we find that the spin-spin correlation 



In 



(g-l)A 



In 



(30) 



1 



- 1 



coth ■ 



Jr 



qiq-2)Jc 



■r 



(31) 



Thus, we find that the correlation length ^ increases as the critical point is approached 
according to ^ ~ t~'^, with critical exponent u = 1. It is interesting to note that the 
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correlation length ^ shows interesting scaling and singular behavior near the critical point. 
While the Ising model on the Bethe lattice exhibits in general mean-field like phase transition 
with "classical" exponents, the critical behavior of the correlation length near the transition 
point coincides with the correlation length behavior in one- dimensional chain with critical 
exponent u = 1, which differs from its "classical value" u = 1/2 |T^. The similar behavior 
of the localization length associated with density-density correlator can be observed in the 
Anderson model on the Bethe lattice IHl. 



If h and t are both sufficiently small in the critical region, then based on Eqs. (122|), (|28 



and ( PID the general behavior of the correlation length ^ should be described by a scaling 
function F 



(32) 



where F{x) = (^fi + /; 



2 



fl 



with 
qiq - 2) 



In 



2(g_l)-g-2 

For a small h > and t ^ 0, we obtain from Eq. 



and /I = 9 



qiq -2) 



(q-iy 

the result: 



(33) 



i.e. the critical exponent is 2/3. 

The bulk susceptibility per lattice site x or the linear response against the field is derived 
from Eq. ( ^1]) as x = dmi/dh 



X 



dm I 
~dh 



Xo 



1 - 



(g- 1) (2 cosh (2 J) -x 
2 sinh (2 J) 



X 



1 -1 



Xo 



l-(g-l)A' 

where xo is nonsingular part of the magnetic susceptibility and given by 



(34) 



Xo 



2e 



_2j(2cosh(2J) - X - x-^){2exp{2J) -x-x 
sinh (2J)(x + a;-i -2exp(-2J))2 



= (l-m?)(l + A). 

Thus, the magnetic susceptibility x can be written finally in the simple form 

(l-ml){l + X) 



X 



l-(g-l)A 



(35) 



(36) 



By means of the fiuctuation relation, x = — ml), we can recover Eq. (|36D 

through Eq. (|27|). To proved this statement let us consider in more detail the fiuctuation 
relation 



^ " N^oo ^ ^2 (< > - <Si>< Sj >), 



(37) 
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where the sum goes over all pairs of sites on the Bethe lattice and Ns is the total number of 
sites. To carry out the summations in Eq. (37), we first note that by definition all sites on 
the Bethe lattice are equivalent and consequently 

< Si > = < So >, 

< SiSj >= Ns^^K SoSj > . 
ij j 

Using these relations we may rewrite Eq. (^7p as 



X 



j 



m 



where cj is the number of sites which is j steps away from the central site 



c, = q{q - iy-\ 



(Co = 1), 



(38) 



(39) 



r(j) is the spin - spin correlation function given by Eq. (]27|) and mi =< sq > is the 
magnetization of the spin in the central site. 

Substituting Eqs. d^) and into Eq. (H), we find 



X = jim (1 - ml) [1 + Y: q{q - iy~'\'] 



lim (1 — mf] 



1 + A 



l-(g-l)A 



-q\ 



n-2 ■ 



(g-l)A 



(40) 



It is now clear that the susceptibility x diverges for {q — 1)X > 1 and given by Eq. (pGf ) for 
(g — 1)A < 1. iFrom Eq. ( pUf ) we can easily established the relation between susceptibility 
X and the correlation length ^ in the critical region 

IV. SPIN - 1 ISING MODEL 

Let us consider, for example, a spin - 1 Ising model, which is known as the Blume- 
Emery-Griffits (BEG) model |31]. The BEG model on the Bethe lattice was studied in Refs. 
]. The model has played an important role in the development of the theory of tricritical 



phenomena 33 



The Hamiltonian of the spin-1 Ising model on the Bethe lattice is given by 
-PH=JY1 SiS, + Kj2 sy^-AY^s^.+hY^Si, 



(41) 



<ij> 



<ij> 



where (3 = {kBT)~^ and Sj = +1, 0, —1 is the spin variable at site i. The first term describes 
the ferromagnetic coupling ( J) between the spin at sites i and j and the second term describes 
the biquadratic coupling {K). Both interactions are restricted to the q nearest neighbor pairs 



10 



of spins. The third term describes the single ion anisotropy A and the last term represents 
the effects of an external magnetic field (h). 

This model has two order parameters, one is the thermal average of the total spin rrii =< 
Sq > and the other is the quadrupolar moment m2 =< Sg > which reflects the possibility of 
phase separation. These order parameters are expressed by 



nil 



exp {h — A)y'^ — exp {—h — A)x'^ 
1 + exp {h - A)yi + exp {-h - A)xi 



(42) 
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exp {h — A)y'^ + exp {—h — A)x^ 
1 + exp {h - A)yi + exp {-h - A)xi ' 



(43) 



iFrom Eq. (|T^) we can obtain the following expression for the first - neighbor spin- spin 
correlation function < SqSi > 



< SqSi >-- 



2x'^y'^e 



-yJ-J+K~2A) 



1 + 2(e^?/^ + e-'^a;^)e-'^ + {e'^'^y^^ + e-2/»a;27)e('^+^-2A) _ 2x^y^e(--^+^-^^) 

(44) 



where 



X = lim 



9n[ 



Let us introduce the new variables 



and y = lim ~^ 



x-y 
V = — - — and u 



X + y — 2 



then we obtain 



mi = —V 



u{l + a) + a 



(45) 



and 



m2 



< SqSi >-- 



u{u + 1) + aw 2 
6 + ^2 + at>2 

_ hu^ + a^{h + l)v'^ 
ah{h + + aw 2) 



(46) 



(47) 



where m, v are the solution of the recursion relation given by Eq. (|T3|) in the thermodynamic 
hmit (A^ oo) 



exp {2h) 



u — av f u + 1 + V 



u + av \u+l — V 



q-l 



(48) 



exp (2A) 



4(6 -m)2 
■^2 — a2t>2 



(m + 1)2 



(49) 
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and b, a are the following constants 

b = exp (K) cosh (J) — 1, 

exp (K) cosh (J) — 1 
exp (K) sinh (J) 

It follows from Eq. (^) that the spin - spin correlation function for spin - 1 Ising model 
can be written as 



m=mi + Y:A,^-^j . (50) 

where Ai and A2 take the form (see Appendix) 

^ _ A3(m2 - ml) - Ai(< spSi > -ml) 
A3 — A2 

and 

^ _ A2(m2 - m^) - Ai(< sqSi > -mf) 
A2 — A3 

with mi, m2 and < sqSi > given by Eqs. (^Sf), (^6]) and (|47|), respectively. The eigenvalues 
Ai, A2 A3 of the symmetric 3x3 matrix V 

/ Vn V,2 Vi3 

V = V,2 V22 V23 
\ Vi3 V23 V33 

can be obtained from the characteristic equation 

(V^ll - X){V22 - A)(V33 - A) + 2VuV,3V23 

= Vi,{Vu - A) + V,liV22 - A) + V.liV^s - A). 
The elements of the matrix V are given by 

Vss' = exp (^Jss' + Ksh'^ - A^^^^ + [x{s)x{s')]'^ , (51) 

where s and s' may independently take values +1, 0, —1. 
Using Eqs.(|D, (H) and (|l|) we find that 

Ai = ^, X2,3 = {C±VC^^) Xi, (52) 

where 

^ u ^ [ab + a + b){u'^ + u — av'^) 



2b 2ab[{u + 1)2 - v'^ 

and 

(6 — m)(m^ — a^v^) 
" afe[(M + 1)2 - t;2] ■ 

12 



The global phase diagram of spin-1 Ising model on the Bethe lattice has been studied 
in detail in the Refs. The A - line of the phase transition in the (J, K and A) space 

given by following conditions: 



/i = and exp (A/ 



^{h-Uc) 



Ur 



1) 



(53) 



where 



Ur 



In terms of the T, Aa/J and K/J, Eq. (^) of the A - line implies a relation T = 
Tc{A\/ J, K / J) which locates the critical temperature as a function of Aa/J and K/J. 

The critical line starts at Aa — oo, Tc/J = |ln^^, which correspond the critical 
temperature of the spin - 1/2 Ising model. We note that for A — oo, the state Sj = is 
suppressed, and the Hamiltonian (^) reduces to the spin - 1/2 Ising model with interaction 
J and external magnetic field h. For certain values of K/ J the system possesses a tricritical 
point at which the phase transition changes from the second to the first order. A tricritical 
point satisfies the following equation 



In particular 



and 



Uc + i o , g - 3 1 

P - Uc 2q Uc 



1 q{3q - 2) K 
-in — ; — , when — 

2 3(g-2)' J 



(54) 



1, -g + 6 + V49g^ - 36g + 36 , K 

- In ; , when — 

2 6(g-2) ' J 



3. 



Let us now consider the general behavior of the spin-spin correlation function in the 
critical region. First consider the case h = {v = 0), then we have 



T(n) = + 



u^ + b 



with 



u 



Ai a{u + 1) 

On the critical line (T = Tc, h = and A = Aa) given by Eq. (|53|) , the spin-spin 
correlation function can be expressed as 



r(n) 



Uc{Uc + 1) 



Ur 



aiuc + 1 J 



a{q-l) 



a? + b{q — 1 — a 
By analogy with spin - 1/2 case we can write 



-(g-1)"". 



(55) 
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r{n) 



q(9 - 1) 
+ b{q — 1 — a 



-n 



(56) 



where d is dimension of the Bethe lattice. For T^/ J = | In ? obtain T{n) = (g — 1) 
which corresponds the case of the spin-1/2 Ising modeL In the critical region, T{n) has an 
asymptotic decay of the form 



Tin) 



exp (-|) 



(57) 



where is correlation length and is given by 



In 



1 



n -1 



(g-l)(A2/Ai) 



Near the critical A - line, setting u — Uc = {uc + 1)6, v = {uc + 1)^ and T = Tc(l + t), 
the h, A and A2 may be expanded, for small 6, e and t as 



h = 'j£[hQt — (a — 1)6 + 



(7^ - 1) 



+ (a^ - 1)52 - (fa - l)^^^ + 



y-1 



(58) 



A - Aa = -Aot + (7 - a - 6)5 + + 



7-1 2 , - 6^ - 7 



5^ — (7a — 1)75^(5 + 7 



7^-1 



(59) 



where 7 = g - l,a = (m^ + l)/uc, 6 = (u^ + l)/(/3 - Uc), 

Jc(exp ~ cosh Jc) — Kc sinh 



/in 



and 



(exp ii'c cosh Jc — I) sinh 

Jc sinh Jc + i^c cosh Jc 

exp Ac 



exp Kc cosh Jc — 1 

Consider the case, when h = and A = Aa- /,From Eqs. (|58|) and ( ^9]) we obtain 

Cit = C26 + 0^6"^ 



(60) 



with 



and 



ci = Ao + 



37 , . Uc+l ^1 7-2 

/^O, C2 = 7 - 1 - 7^ h 



2(7+1) 



C3 



a2 - 62 _ ^ 



/5-Mc Mc 2(7 + 1) 

37 1 37(72-7-3) 1 



2 7 + lMc 4(72 -1)(7 + 1)^2 ■ 

It is easy to see from Eqs. (plf ) and (|60D that in all points on the A-line t ~ 5, except for 
the tricritical point, where t ^ 6"^. 

Using Taylor expansion of the expression for the correlation length ^ by small e, 6 and t 
and Eqs. ( p^ and (§3), we find that the correlation length ^ increases as the critical point 
is approached according to ^ ~ i^^, with critical exponent z/ = 1, everywhere on the A-line 
except the tricritical point, where ^ ~ t^^/2 ^ith tricritical exponent Ut = 1/2. 
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V. SUMMARY AND DISCUSSION 

Let us now briefly summarize our results. In this paper we consider the most general 
spin - S model on the Bethe lattice in the external magnetic field and use the transfer - 
matrix approach to derive the spin-spin correlation function. The general spin - S model 
indues the spin-1/2 simple Ising model and BEG model as special cases. /^From the exact 
spin - spin correlation functions T{n) of the Bethe lattic Ising model and BEG model in an 
arbitrary magnetic field h and temperature T, the correlation length ^ has been determined 
analytically. In the critical region the correlation length ^ of the simple Ising model is 
inversely proportional to the distance from the critical point. Such singular behavior 
coincides with the correlation length behavior in the one - dimensional chain. We obtain also 
that near the transition point the magnetic susceptibility is proportional to the correlation 
length ^. 

Recently, Gujrati [ll|] showed that in many cases the behavior on Bethe or Bethe-like 



lattices are qualitatively correct even when conventional mean-field theories fail. By a proper 
choice of these lattices, it is possible to satisfy frustrations, gauge symmetries, etc., which 
are usually lost in conventional mean-field calculations, because of the lack of correlations. 
Such correlations are present on the Bethe-like lattices, and in this paper we have given the 
exact expression for such correlations. 

It should be noted that we can obtain the proper singular behavior of ^ because we 
have used the proper thermodynamic limit (A^ —>■ oo) to obtain the recursion equations 
and correlation function, i.e. Eqs.(|T3p and (|Hp for the most general spin - S model on the 
Bethe lattice and we have demonstrated the crucial role of Bethe lattice dimensionality in 
determining critical behavior of the correlation length. 

In conclusion, it must be remarked that the transfer - matrix methods discussed in 
this paper can be extended without difficulty to obtain correlation functions with singular 
correlation length for other spin models, e.g. Potts model, multilayer Ising model, Ising 
model with competing nearest-neighbor and next nearest-neighbor interactions, etc, on the 
Bethe and Bethe-like structures. This approach should be applicable for gauge models on 
generalized multi-plaquette hierarchical structures as well. 
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APPENDIX 

The coefficient Ak {k = 1, 2, 25*) can be obtained by solving the following system of linear 
algebraic equations 



f g^ = Ai + A2 + ... + A2s 

g2 = Aih + A2I2 + ••• + ^25^25 

[g2s = A,lf-' + A2lf-' + ... + A2slfs-\ 
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where 

Qk =< SoSk-i > -ml ior k^l,2,...,2S 

and Ik = Afc+i/Ai. 

Let us introduce the elementary Lagrange interpolation polynomials 



^i(0 = n rzf = EM'"\ ^ = 1,2, ...,25, 

37^» 



which satisfy Li{lj) — Sij, where Sij is the symbol Kronecer. Now, it is evident that 

25 

^k^Yj^kiQi, A; = 1,2, ...,25". 

The elements of a^j are 

an, = 

where is 25 — 1 degree polynomials defined by 

25 

Pk{l) = \[il-h), k^l,2,...,2S 

i=l 

and Fj{k) is the elementary symmetric function in 25 — 1 variables li,l2, h-i, h+i, 

F2S = 1, 

-^25-1 = ^1 + ■■■ + h-1 + + ■■■ + ^25, 

Fi = h-.-lk-ih+i-'-hs- 
Thus, the coefficients will take the form 



riifii^fc — k 

Examples: 
1. S=l/2 



■A-i = gi =< si > —ml = 1 — ml- 



For spin - 1/2 Ising model < Sq >= 1. 
2. S=l 



h-h 

\^{m2 - ml) - Ai(< sqSi > -ml) 
A3 — A2 



and 
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A.2 ■ 



A2(m2 - m\) - Ai(< sqSi > -ml) 
Ao — As 



where m2 =< Sg >. 
3. S=3/2 



with 



< soSn > -ml = Ail'l + A2I2 + Asl^, 
khgi - {k + h)g2 + 5'3 



{k-li){h-k) 



and 



A2^ A^ih^h), A^^ A^{h^h), 

where gi =< Sq > —ml, g2 —< sqSi > —ml and gs —< S0S2 > -ml- 
4. S=2 

< soSn > -ml = A^q + A2q + Aa/^ + AJ2- 
Figi - F2g2 + -Fsfl'a - F^g^ 



with 
and 



{l2-h){h-h){h-k) 

Fi = hhh, F2 = I2I3 + hh + hh, Fs = I2 + h + h, F4 = 1 
Ak = Aiih^k), ior A; = 2,3,4, 



(63) 



where gi —< Sq > 92 —< sqSi > —ml, gs —< S0S2 > —ml and g^ —< so^s > 



-ml- 
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